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Abstract. We give an algorithm for efficient step size eontrol in numerical integration of non- 
stiff initial value problems, based on a formula tailormade to methods where the numerical solution 
is compared with a solution of lower order. 

Key words. Ordinary differential equations; Step size selection; Efficiency 



1. Introduction. We discuss here step size control in numerical integration of 
ordinary differential equations, for the purpose of optimizing performance in terms 
of precision and computer time. Thus the aim is not to estimate or control absolute 
errors, although upper limits can be given. Only non-stiff initial value problems are 
considered. 

A new formula for efficient step size control has recently been proposed . This 
formula has the desirable property that the value it gives for the next step size 
is independent of the present step size /i„ in the asymptotic limit /i„ — ^ 0. For a 
one-step, variable step size method of order {p,p—l) that in addition to the primary 
integration method of order p, also uses a secondary method of order p — I, the step 
size control algorithm using this formula can be written: 



(1-1) --/t-r-L '-'t III; 

\e\h„ 

(1.2) L : Xihn < hn+i < A2/i„ 

(1.3) A : \e\K < <Jt 

Here e is the difference between the two numerical solutions in one single step of length 
hn and r is a constant tolerance. We limit the increase and decrease of the step size by 
the limitation criterion L, and reject steps that do not satisfy the acceptance criterion 
A. a, Ai and A2 are parameters. We show that the above algorithm is more efficient 
than the standard algorithm known from the literature B: 



(1.4) 
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(1.5) L : Xihn < hn+l < X2hn 

(1.6) A : |e| < ar 



One way to measure efficiency for one-step methods is to determine the number of 
function evaluations for a given global error. We use the software package DETEST 
, applying the embedded Runge-Kutta method of Dormand and Prince Q , of order 
{p,q) = (5,4). First we find the values of a, Ai and A2 that give the least number 
of function evaluations for a global error of 10^^, 10"^, 10^^ and 10^^. This is 
done simultaneously for 10 of the problems (group I) in DETEST using a Levenberg- 
Marquardt method with a mixed quadratic and cubic line search procedure We 
then compare the two algorithms with this set of parameters on 10 other problems 
(group II) from DETEST. The computer used was a Cray Origin 2000. 
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2. Numerical Results. Group I is defined arbitrarily to be the problems Al, 
A3, A5, B2, B4, CI, C3, C5, E2 and E4, and group II to be the problems A2, A4, Bl, 
B3, B5, C2, C4, El, E3 and E5, in DETEST §. We omitted the problems in class 
D and F, since DETEST only gave results for a limited set of tolerances in these (the 
tolerance is varied automatically by the program). The parameters that resulted in 
the least number of function evaluations for group I were: 



(2.1) 
(2.2) 



(T = 6.70 Al = 0.67 A2 = 5.00 (formula 1.1), 
cr = 5.50 Ai = 0.26 A2 = 4.00 (formula 1.4). 



The new formula gave a lower number of function evaluations in 50.9% of the cases in 
problem group I (for which cr, Ai and A2 were optimized), and in 64.4% of the cases in 
problem group II. Only calculations giving global errors in the range (10^"^, 10^^) were 
count ed, s ince DETEST did not give results outside this range for some problems. In 
Table 2.1 we have displayed the mean ratio of the number of function evaluations of 
the new formula to the number of function evaluations of the standard formula. The 



standard formula was tested with the optimized set of parameters (equation 2.1) and 
the recommended set of parameters: 



(2.3) 



(7=1.20 Ai=G.50 A2 = 2.GG, 



Table 2.1 

The mean ratio of the number of function evaluations of the new formula to the number of 
function evaluations of the standard formula. 



Expected 
global 
error 



with equation and 2.2 
Group I Group II 



Mea ri ra tio of function evaluations 



with equation 2.1 and 2.3 


Group I 


Group II 


0.8794 


0.8747 


0.9222 


0.9011 


0.9281 


0.9266 


0.9370 


0.8797 


0.9397 


0.8705 


0.9434 


0.8814 



10 
10 
10 
10 
10 
10 



0.9575 
1.0152 
1.0166 
1.0014 
1.0000 
0.9947 



1.0103 
1.0448 
0.9696 
0.9615 
0.9518 
0.9545 



3. Discussion. A new formula for step size selection in numerical integration 
of non-stiff initial value problems has been tested on 20 initial value problems. It is 
found that this formula on the average is more efficient than the standard step size 
selection formula. I thank the authors of DETEST. 
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